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We discuss the impact of the true value of the CP phase on the mass hierarchy, CP violation, and 
CP precision measurements at neutrino superbeams and related experiments. We we use a complete 
statistical experiment simulation including spectral information, systematics, correlations, and de- 
generacies to produce the results. However, since it is very complicated to understand the results 
in terms of a complete experiment simulation, we show the corresponding bi-rate graphs as useful 
tools to investigate the CP phase- dependencies qualitatively. Unlike bi-probability graphs, which 
are based upon oscillation probabilities, bi-rate graphs use the total event rates of two measurements 
simultaneously as a function of the CP phase. Since they allow error bars for direct quantitative 
estimates, they can be used for a direct comparison with a complete statistical experiment simula- 
tion. We find that one can describe the CP phase dependencies of the mentioned measurements at 
neutrino superbeam setups, as well as one can understand the role of the sg^AmJ^-degeneracy. As 
one of the most interesting results, we discuss the dependence of the CP precision measurement as 
a function of the CP phase itself, which leads to "CP patterns" . It turns out that this dependence 
is rather strong, which means that one has to be careful when one is comparing the CP precisions 
of different experiments. 



PACS numbers: 14.60.Pq 



I. INTRODUCTION 



The leading atmospheric and solar neutrino oscilla- 
tion parameters are now measured with unprecedent- 
edly high precisions. The Super-Kamiokande 0, H 0, 
SNO |j, la, 13 j and KamLAND Q experiments have es- 
pecially contributed to these achievements. Furthermore, 
the establishment of the LMA (Large Mixing Angle) so- 
lution makes sure that three-flavor effects are not en- 
tirely suppressed by the mass hierarchy. This means 
that they could be detectable by the next-generation 
experiments - provided that the leptonic mixing angle 
6*i3 is large enough. So far, #13 has been constrained to 
sin 2 20i3 < 0.1 by the CHOOZ experiment [10. In or- 
der to determine #13, the neutrino mass hierarchy, and 
the leptonic CP phase Sqp, conventional beam experi- 
ments [lOt llll | l2l are op erated or built, and superbeam 
experiments El HH El are proposed. As an alterna- 
tive, new reactor experiments with near and far detectors 
could provide some information on #13 on a relatively 
short timescale El 13 El E! After the superbeam 
era, future neutrino factories could find #13 > and test 
leptonic CP violation even below sin 2 29\3 ~ 10~ 3 (for a 
summary, see Ref. ppf). 

The next generation of experiments which are sensi- 
tive to all of the unknown oscillation parameters are 
the neutrino superbeams. In this study, we especially 
consider the JHF (J-PARC) to Super-Kamiokande (or 
Hyper-Kamiokande) El and NuMI El 

off-axis super- 
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beams, since specific LOIs or proposals exist for these 
experiments. Unfortunately, the potential to access 
#13, the neutrino mass hierarchy, and the leptonic CP 
phase ^cp simultaneously implies that superbeams are 
facing a high-dimensional and very complex parameter 
space. In particular, the measurements are spoilt by 
multi-parameter correlations l2ll and intrinsic degenera- 
cies, which are the (<5cp,#i3) u3- sgn(Am| 1 ) [23], and 
(023,tt/2 — 823) de gene racies leading to an overall 
"eight-fold" degeneracy |25f. In principle, one could ei- 
ther combine two first-generation superbeams, or one su- 
perbeam with a reactor experiment to resolve these de- 
generacies for sin 2 26 13 > 10" 2 [H El [H [H [H H S 
l31| . since these experiment types are proposed for com- 
parable timescales and sin 2 2#i3-values. If one wants to 
resolve the degeneracies with superbeams in the range 
sin 2 2#i3 ~ 1CU 2 — 1CU 3 , superbeam upgrades, such as 
the JHF (J-PARC) t o Hy per-Kamiokande setup, will be 
required [33, Hj, [34], y^, For a systematic discus- 

sion about possibilities to resolve degeneracies on differ- 
ent scales of sin 2 26>i3, see Ref. j^l and references therein. 

The simulation techniques for future long-baseline ex- 
periments are very advanced now. We use in this study 
a complete statistical simulation including spectral infor- 
mation, a realistic detector simulation, systematics, cor- 
relations, and degeneracies. However, the results from 
such a simulation are not always trivial to understand, 
especially because of the impact of correlations and de- 
generacies. Therefore, it becomes more and more impor- 
tant to have tools to interpret these results. Note that 
these tools will never be able to replace the complete nu- 
merical simulations, but they can help to understand at 
least some of the aspects of the results. In this study, 
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we use one of these tools, the bi-rate graphs, which are 
useful to understand the CP phase-dependencies of many 
measurements. The impact of Sep is rather large in the 
considered experiments, since many of the problems with 
correlations and degeneracies originate in the correlation 
with 8qp- We will mainly discuss the experiment perfor- 
mances as function of the "true value of Sep" ■ We refer 
to the true value of Sep as the simulated value, which is 
realized by nature and which is, up to know, completely 
unknown. 

Bi-rate graphs are based upon the so-called "bi- 
probability" graphs used in Refs. |H H IS |U HI EJ 
|42| . In these graphs, the CP phase-dependence results 
in ellipses in the space of the oscillation probabilities 
of the two considered oscillation channels. Comparing 
the ellipses for different values of sin 2 28i 3 and for differ- 
ent degenerate solutions leads to qualitative descriptions 
of many CP phase-dependent effects. In this work, we 
use the bi-rate space of the total charged-current event 
rates instead of the oscillation probabilities, such as it 
was first suggested in Ref. |23j . Since error bars can 
be used in this approach, semi-quantitative statements 
can be directly compared to a complete statistical anal- 
ysis including systematics, multi-parameter correlations, 
and degeneracies. Because the off-axis superbeams are 
"quasi-monochromatic" beams due to their narrow beam 
spectrum ( "narrow band beams" ) , spectral information 
plays only a minor role. Therefore, we will demonstrate 
that one can describe many properties of the results of 
complete statistical superbeam simulations with the help 
of bi-rate graphs, and thus the underlying bi-probability 
picture. 

Besides the combination of superbeams, we will dis- 
cuss the results including a new proposed reactor experi- 
ment in some cases, because we will observe that the 
results have a similar structure. We do not investigate 
future neutrino factories or other high-precision instru- 
ments, because the interpretation of the bi-rate approach 
is much more complicated for beams with good spectral 
information. Note that we treat the chosen experiments 
as examples in this study, which means that we do not 
want to give selection criteria for specific setups. The 
purpose of this study is to understand the CP phase- 
dependencies of complete statistical simulations in terms 
of bi-rate graphs phcnomcnologically. In fact, we will ob- 
serve that the qualitative main results to not depend on 
the chosen setups. 

This study is organized as follows: In Sec. [HI we an- 
alytically introduce the bi-probability picture and sum- 
marize some of its basic properties. This section can be 
skipped by the reader familiar with the bi-probability ap- 
proach. The complete statistical and bi-rate experiment 
simulations are described in Sec. II I II where the experi- 
ment setups used as examples are introduced, too. As 
the main part of this work, in Sees. IIVI Ivl and IVII the 
CP phase-dependent properties of the mass hierarchy, CP 
violation, and CP precision measurements are discussed. 
Finally, we summarize our observations in Sec. IVIII 



II. GENERAL PROPERTIES OF 
BI-PROBABILITY GRAPHS 

Bi-probability graphs, as they were introduced in 
Refs. |H H3, HEM US El E3, show the probabili- 
ties of two measurements Pi and P 2 simultaneously as 
a function of Sep- For a given set of parameter values, 
the oscillation probabilities Pi and P 2 lead to a point in 
the bi-probability space (Pi,P 2 ). Varying the CP phase 
from to 2tt, one obtains ellipses with certain character- 
istics, which allow to study many properties depending 
on Sep- A very nice introduction to bi-probability graphs 
and their properties illustrated by simple figures can, for 
example, be found in Ref. Using total event rates 

instead of oscillation probabilities, one can also draw bi- 
rate graphs with similar properties. The advantage of bi- 
rate graphs is the possibility to use error bars in order to 
obtain direct quantitative estimates of many effects [2^| . 
Many of the features introduced below can, for example, 
be found in Fig. [21 (left), which may serve for illustration. 
Typical cases for the measurements described by Pi and 
P 2 are: 

1. The probabilities (P^J , PjS) of two different su- 
perbeam experiments, i.e., for two different values 
oiL/E. 

2. The CP-conjugated channels (P Me ,P^ g ) of one su- 
perbeam experiment. 

3. The T-conjugated channels (P^Pe/J for the in- 
vestigation of T violation. 

In this study, we will not consider case 3, since the frac- 
tion of u e produced by pion or kaon decays is rather small 
for superbeams. 

In order to obtain some qualitative analytical under- 
standing of bi-probability graphs, we can write the os- 
cillation probabilities for a small #13 and a small mass 
hierarchy parameter a = lS.m\ 1 /\/S.rn\ 1 \ as [30l l39| 



Pi=X. 
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In these equations, i refers to the considered measure- 
ment i € {1,2}, Li to the baseline, Ei to the neutrino 
energy, and X* to the electron number density in mat- 
ter. In addition, s a & = sm6 a b and c a b = cos 9 a b- The up- 
per signs in the individual terms correspond to neutrino 
running and a normal mass hierarchy. 1 For the inverted 
mass hierarchy, all the double signs change to the lower 
ones. For antineutrino running, the double signs in the 
\A\ 3 =F cti-Li | -terms and in Eq. Q change to the lower 
ones. Note that using the inverted mass hierarchy as 
well as antineutrino running causes double sign changes 
at the respective places, i.e., some signs are changed back 
to the original ones. 

In general, Eq. (JIJ has the same structure as the orig- 
inal expansions in Refs. [iM E3, l45| : the X^-term is pro- 
portional to sin 2 2# 13 , the Y^-term to sin2# 13 a, and the 
P®-term to a 2 . Thus, the CP effects are suppressed by 
the small #13 and the mass hierarchy parameter a. One 
can also proove that the CP trajectories are elliptic and 
that the minor and major axes of the ellipses are in case 2 
in vacuum always parallel to the diagonal and its orthog- 
onal for equal L/E (c/., Appendix of Ref. 23]). However, 
for different values of L/E, the pre- factors (offsets) in Pi 
and Pi can be different [c/., Eq. and the symmetry 
axis is not the diagonal anymore. A similar effect occurs 
for bi-rate graphs, because cross-sections, efficiencies etc. 
lead to different pre- factors even for the same L/E. 

From Eq. QJ, we can derive some basic properties of 
bi-probability graphs. For fixed oscillation parameters 
(except from <5cp)j this set of two equations gives the 
parameterized bi-probability trajectory in (P\, Pi)-sj>ace 
as a function of Sep. It is obvious that the trajectory 
becomes degenerate if Pi and P2 have the same phases 
in the cosine term in Eq. In this case, the ellipse 
simply collapses to a line. In order to create a phase 
shift between P\ and Pi in this term, one can, for exam- 
ple, use a different L/E for Pi and Pi (different exper- 
iments, case 1), or neutrinos in P\ and antineutrinos in 
Pi (CP-conjugated channels, case 2). The special case of 
a phase shift 7r/2 results in a circle, such as it is obtained 
for one experiment in the oscillation maximum and the 
other in the oscillation minimum in case 1. The positions 
of the ellipses are determined by the offset proportional 
to sin 2 (except from the small contribution from P®), 
and the lengths of the axes of the ellipses are proportional 
to sin ^13. Thus, it is usually quite illustrative to show 
many different ellipses for different values of sin $13 , such 
as it is done in Fig. [21 left. In addition, the solutions for 
the normal and inverted mass hierarchies are normally 
drawn in each figure in order to discuss the effects of 
the sgn(Am§i)-degeneracy j^. Though the ellipses for 
the normal and inverted mass hierarchies are identical 
in vacuum, they move away from the diagonal (for equal 



Note that we only consider the oscillation channel Vj_ L — » v e in 
this work. 



L/E's in Pi and Pi) in opposite directions for an increas- 
ing matter potential of one experiment. This asymmetry 
comes from the different sizes of the |A^ 3 =F o-iLi |-terms 
for the different mass hierarchies and experiments, which 
modify the offsets [c/., Eq. J5J)]. Since one point in the bi- 
probability space corresponds to several sets (#i3,<5cp), 
i.e., different ellipses, it makes sense to show the area 
of all possible bi-probability trajectories generated by all 
possible values of #13 for a specific mass hierarchy, which 
is often called a "pencil" (c/., Fig.|3 left). Thus, a unique 
determination of the set (#13, <5cp) is sometimes impossi- 
ble because of the correlation between #13 and <5cp- In 
addition, if the two pencils for the positive and negative 
mass hierarchies overlap, then the situation becomes even 
more complicated because of the opposite-sign solutions. 

In Eq. Q), the phase 4>i = A^ 3 /2 in the cosine term 
is 7r/2 at the first oscillation maximum. Thus, close to 
the first oscillation maximum, we can make the approx- 
imation cos(<5cp ± 4>i) — cos((5cp ± 7r /2) = Tsin^cp- 
Assuming that Pi and Pi represent case 1 (two different 
experiments both using neutrinos), we can make some 
simple observations, which will be useful for the follow- 
ing sections: 

Change of L/E: If <\>\ = <fii, i.e., for the same L/E, 
the ellipses collapse to lines. Moving from <f>i < (f>i 
to <f>\ > 4>i (or vice versa) changes the orientation 
of the CP trajectories, as it can be easily seen from 
Eq. 0. 

CP phase closest to origin: For cf>i ~ tfii ~ ir/2, i.e., 
both experiments operated close to the oscillation 
maximum, cos(<5cp + tt/2) is — sin<5cp for the nor- 
mal mass hierarchy. Thus, the CP trajectory is 
closest to the origin for Sep — n/2 and farest away 
from the original for Sqp — 3ir/2. For the inverted 
mass hierarchy, the additional sign change in <f>i is 
compensated by the sign change in Yi [c/., Eq. (J3J], 
which means that we have the same behavior. 

Inverted mass hierarchy: Since cos(£cp — (f) — 
cos(— <5cp + 4>) and the closest and farest values 
(from the origin) of the CP trajectory do not 
change for the inverted mass hierarchy, the CP tra- 
jectory changes its orientation. In addition, the 
sin (|A l 13 =F a i L l \/2) / (\A\ 3 =p a,ij|/2) -terms in A; 
and Yi [c/., Eqs.|21and|3| become somewhat smaller 
for the inverted mass hierarchy close to the oscilla- 
tion maximum, which means that the ellipses move 
closer to the origin and become somewhat smaller - 
provided that matter effects are present. The same 
holds for the Pf-ter m m Eq. (£). Matter effects 
also cause the asymmetry between the positive- and 
negative-sign pencils coming from different A^'s in 
Eq. ©. 

Effect of the solar Am^: The size of the solar Am^ 
changes the offsets [c/., Eq. Q}] as well as the axes 
lengths [c/., Eq. (J3J] of the ellipses. Thus, the larger 
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the value of Ato^i is, the more will the ellipses be 
extended. 

Very small values of #13: For very small values of #13, 
the smallness of #13 dominates the splitting be- 
tween the normal and inverted hierarchies, since 
in this limit the matter effect becomes a second or- 
der effect. Thus, the ellipses for the normal and 
inverted mass hierarchies strongly overlap. 

Similar observations can be derived cases 2 and 3, which 
we however will not discuss further. However, it is in- 
teresting to note that for case 2 the CP conserving solu- 
tions and 7r lie on the line P\ = P% in vacuum because 
of the identical oscillation probabilities for neutrinos and 
antineutrinos. In addition, in this case usually an equal 
L/E (same experimental setup) for Pj and Pi is assumed. 

Since Eq. Q includes three terms, which can be of sim- 
ilar order of magnitude, other correlations than between 
<5cp and #13 may be present and influence the measure- 
ment. So far, we have assumed that all parameters except 
from <5cp and #13 are fixed. As long as the correlation be- 
tween Sep and #13 dominates, the bi-probability picture 
will give the proper results. However, if there is a sub- 
stantial contribution of another correlation, then there 
may be substantial differences to a complete statistical 
analysis including multi-parameter correlations. For ex- 
ample, the size of Amli, and thus P®, is obtained from 
solar neutrino experiments. Since it will not be exactly 
known at the time of the superbeam operation, the offset 
of the ellipses may vary within certain limits and lead to 
another error source. This example illustrates why we 
will find limitations of the bi-probability formalism in 
some cases. 



III. BI-RATE PICTURE AND COMPLETE 
EXPERIMENT SIMULATION 

In this section, we first of all introduce the bi-rate 
approach as modification of the bi-probability picture, 
and discuss its limitations. Then, we introduce the com- 
plete statistical simulation techniques shortly, and point 
to earlier studies where they are described in greater de- 
tail. Eventually, we give the details of the experiment 
examples used in this study, as well as the oscillation 
parameters used. 

In order to directly compare the bi-probability graphs 
to the results of a complete statistical analysis, we use 
the bi-rate picture. Instead of the bi-probability space, 
it shows the ellipses in the bi-rate space of the total (inte- 
grated) charged-current event rates. The main advantage 
of this approach is that statistical error bars can be used 
as a direct measure for the statistical significance of a 
specific result. Because of the off- axis technology to pro- 
duce narrow-band beams, the superbeams usually con- 
tain only very little spectral information. Thus, the total 
rates of the considered experiments carry the dominant 



physics information. Especially, the charged-current ap- 
pearance rates dominate the measurements of the most 
interesting parameters sin 2 29 13, the mass hierarchy, and 
<5cp • Unfortunately, a small fraction of the detected elec- 
tron neutrino (or antineutrino) signal is not originally 
part of the charged-current appearance rates and is of- 
ten called "background" . For the bi-rate approach, we do 
not use the background events, since they do not contain 
the interesting information (of course, we include them 
in the complete statistical simulation) . Similar to the bi- 
probability graphs, one can show the total event rates in 
bi-rate space for any combination of two measurements. 
An example for a bi-rate graph can be found in Fig. |2 
which we will use later. In this and similar figures, several 
ellipses for selected different values of sin 2 29is are shown, 
as it is described in the figure captions. The smaller 
sin 2 2#i3 is, the smaller are the total rates and the closer 
are these ellipses to the origin. The different values of 
5cp are marked by symbols, where the triangles refer to 
the CP conserving solutions. The total event rates for 
the normal mass hierarchy are shown with solid curves, 
where the regions of all possible ellipses ("positive-sign 
pencils" ) are marked with the dark (green) shading. The 
total rates for the inverted mass hierarchy are shown with 
dashed curves, where the regions of all possible ellipses 
( "negative-sign pencils" ) are marked with the light (yel- 
low) shading. At selected points, the statistical error bars 
for the total rates are shown at the 90% confidence level. 
Note that they are only drawn for the directions along 
the axes in order to obtain an idea of their relative size. 

As we will see later, many properties depending on Sep 
can be understood in terms of bi-rate graphs. However, 
this does not mean that a bi-rate graph can replace a 
complete experimental simulation. First of all, the exper- 
iment might be better than what one can see from the 
bi-rate picture, because there is at least some spectral 
information and in addition some information from the 
disappearance channels. For example, the sgn(Am 2 1 )- 
degenerate solution may actually not appear under the 
chosen confidence level. On the other hand, the experi- 
ment might be worse than the bi-rate interpretation, be- 
cause the latter does not include systematics and other 
correlations than the one between sin 2 2#i3 and Sqp ■ 
However, we will be able to describe many CP phase- 
dependent properties of a complete simulation qualita- 
tively, especially in the regions of the parameter space 
where the total appearance rates carry the dominant in- 
formation. Therefore, the bi-rate picture can be used as 
a tool to understand many basic features of very compli- 
cated experiment simulations. 

In order to demonstrate many effects quantitatively as 
realistic as possible, we use in parallel to the bi-rate ap- 
proach a complete statistical analysis including systemat- 
ics, multi-parameter correlations, and degeneracies. The 
corresponding analysis techniques are introduced and dis- 
cussed in the Appendices of Ref. (2l| (general techniques) 
and Ref. [23 (specifics of the superbeams), as well as in 
Ref. 0| (reactor experiments). In order to use all avail- 
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JHF-SK 


NuMI 


JHF-HK 


Beam 


Target power 


0.77 MW 


0.4 MW 


4MW 


///^-running 


v only 


u only 


v and v 


Running period 


5 years 


5 years 


2 years years (p) 


Detector 


Technology 


Water Cherenkov 


Low-Z calorimeter 


Water Cherenkov 


Fiducial mass 


22.5 kt 


50 kt 


1 000 kt 


Baseline 


295 km 


712 km 


295 km 


Off-axis angle 


2° 


0.72° 


2° 


Overall facts 


Mean L/E 


385 km GeV -1 


320 km GeV -1 


385 km GeV" 1 


Mean energy 


0.76 GeV 


2.22 GeV 


0.76 GeV 


Signal events 


132.8 


334.4 


12 594 0) + 9 239 (P) 


Background events 


22.8 


57.1 


2159 {y) + 3 352 {v) 



TABLE I: The two first-generation superbeams JHF-SK and NuMI, and the superbeam upgrade JHF-HK, as originally given 
in Refs. |T3. IT^j. as well as the standard parameters used in this study. The total signal and background events are calculated 
for the appearance channels and the parameter values sin 2 26^ = 0.1, <5cp = 0, [Am§i| = 2.5 ■ 10~ 3 eV 2 , sin 2 2#23 = 1.0, 
Amii =7.0- 10- 5 eV 2 , sin 20i2 = 0.8, and a normal mass hierarchy. 



able information, we fit the appearance and disappear- 
ance rates simultaneously, where we include the back- 
grounds as described in Refs. 0, |2£j- Therefore, we 
assume that the best measurement of the leading at- 
mospheric parameters will be provided by the experi- 
ment (or experiment combination) itself at the time of 
the analysis. As external input, we assume the product 
Am|x • sin2# 12 to be measured with 15% precision by 
the KamLAND experiment b y t he time the superbeam 
experiments will be analyzed [46L |47| . 

As examples for experiments, we study in this work 
the JHF (J-PARC) to Super-Kamiokande (or Hyper- 
Kamiokande) ^| and the NuMI 0] projects, using the 
beams in the off-axis configurations. We refer to the 
first-generation superbeams with JHF-SK and NuMI 
further on, as well as we use the JHF (J-PARC) to Hyper- 
Kamiokande superbeam upgrade JHF-HK. Their stan- 
dard configurations, as they are proposed in the respec- 
tive LOIs, are summarized in Table [I] Note that we use 
a NuMI detector size of 50 kt (as, for example, in the 
NuMI proposal). We use five years of neutrino running 
time for each of the first-generation beams as well as two 
years of neutrino running and six years of antineutrino 
running for the JHF-HK upgrade in order to have ap- 
proximately equal numbers of neutrino and antineutrino 
events. In addition, we use several modifications com- 
pared to Table HJ In order to be sensitive to the mass 
hierarchy, one needs a longer NuMI baseline |29ll30Ll3^ . 
Therefore, we use a modified NuMI baseline of 890 km 
for the mass hierarchy measurement, which is the longest 
possible baseline for the fixed decay pipe and the con- 
sidered off-axis angle. For the CP measurements, some 
CP-complementary information is needed to resolve the 
correlation between sin 2 26*13 and 5cp ■ This can either 



be a substantial fraction of antineutrino running, or a 
large reactor ex peri ment for a "clean measurement" of 
sin 2 26»i3 EElSEii- Therefore, we either use NuMI 
in the antineutrino running mode only, or the large reac- 
tor experiment Reactor-II from Ref. f° r CP mea- 
surements at the first-generation superbeams. However, 
note that for the first-generation superbeams the detec- 
tion of CP violation at LMA-I will be very hard [29| . 
Nevertheless, we will demonstrate that one could still 
learn something about the CP phase in some cases. 

For the atmospheric oscillation parameters, we use, 
if not otherwise stated, lAmU = 2.5 • 10~ 3 eV 2 and 
sin 2 2023 = 1-0 3]. For the solar oscillation parame- 
ters, we use Am|i = 7 • 10" 5 eV 2 and sin 2 29 12 = 0.8 
(see, for example, Ref. |4g|V Note that the LMA-II so- 
lution is disfavored by the SNO salt-enhanced neutral- 
current data (or, for example, Ref. HH). For 613, 
we only assume sin 2 29\3 < 0.1 9] below the CHOOZ 
bound and for 5cp we do not make and special assump- 
tions, i.e., 5gp € [0, 2?r[. I n ge neral, we include the full 
{Sep, Ota) El, sgn(Am 2 i) and (0 23 , tt/2 - 6 23 ) 
degeneracies, which lead to an overall "eight-fold" de- 
generacy [25j . However, since the energy resolution of 
the superbeams is rather poor and the atmospheric mix- 
ing angle is set to maximal mixing, the (<5cp,^i3)- and 
(023,7r/2 — ^^degeneracies are not relevant for this 
study. 



IV. THE DETERMINATION OF THE MASS 
HIERARCHY 

We define sensitivity to a specific sign of Ato§i (i.e., 
the normal or inverted mass hierarchy) if there is no de- 
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FIG. 1: (Color online) The sensitivity to the mass hierarchy for JHF-SK and NuMI@890 km combined, where both are 
operated with neutrino running only. The symbols on the ellipses correspond to Sep = (filled triangles), <5cp — it/2 (stars), 
Scp = I" (unfilled triangles), and Sep = 3n/2 (diamonds). The ellipses are shown for sin 2 20i3 = 0.01 and 0.08, respectively, 
from the smallest to the largest, and the error bars for the total rates are drawn at the 90% confidence level. The dark shaded 
(green) regions represent the positive-sign pencils with the solid curves, whereas the light shaded (yellow) regions represent the 
negative-sign pencils with the dashed curves. In the right-hand figure, the regions of sensitivity to a normal (sgn(Am| 1 ) > 0) 
and inverted (sgn(Am3!) < 0) mass hierarchy are shown for the complete statistical experiment simulation including all 
information, where the sensitivity is given on the right-hand side of the curves at the 90% confidence level. The horizontal lines 
are marked by the corresponding symbols from the left-hand plot, and the vertical lines correspond to the conservative reach 
in sin 2 2#i3 . 



generate solution with the opposite sign of Am^ fitting 
the simulated rate vector at the chosen confidence level. 
Since especially the sgn(Am§ 1 )-degeneracy [2j| implies 
that changing the sign of Am 2 ^ can be often compen- 
sated by a different value of Sep, it constrains the poten- 
tial of an experiment to determine the mass hierarchy. 
This already indicates that the true value of Sep affects 
the mass hierarchy potential. 

We show in Fig. ^ the sensitivity to the normal and in- 
verted mass hierarchy at the example of JHF-SK and 
NuMI@890 km combined. In this figure, the bi-rate 
graph (left) is compared to the complete statistical sim- 
ulation including spectral information, systematics, cor- 
relations, and degeneracies (right). Note that choosing 
a shorter NuMI baseline would lead to a much weaker 
mass hierarchy sensitivity, as we have discussed in the 
last section. For our choice, the L/E's of both exper- 
iments are very similar, which means that the bi-rate 
ellipses shrink to lines |30( and favor the mass hierarchy 
sensitivity. If the error bars at a chosen point in bi-rate 
(corresponding to a set of true parameter values) touch 
or overlap the opposite sign pencil, then the sensitivity to 
the positive sign of Am^ will be destroyed at the chosen 
confidence level. Note that for a given sin 2 26\3, the true 
value of Scp strongly influences the potential to establish 
the considered hierarchy, since the corresponding point 
may lie closer or farer away from the opposite-sign pen- 



cil. In addition, because of the correlation between Sep 
and sin 2 26x3, any point within the opposite-sign pencil 
will destroy the sensitivity, which does not necessarily 
have the same values of Sep and sin 2 29^. In fact, one 
can identify the CP phase of the fake solution as the CP 
phase of the point closest to the original solution, since 
the total rate differences between those two points are 
smallest. This CP phase of the fake solution does usu- 
ally not change as long as the pencils do not overlap. 
However, if they overlap, then it will start to move (c/., 
Fig. 8 of Ref. 

One can obtain an estimate for the sin 2 2#i3 sensitiv- 
ity reach, i.e., the smallest value of sin 2 29 13 for which 
the considered sign of Am^ can be established, by find- 
ing the ellipse for which the error bars just touch the 
opposite-sign pencil (for a fixed value of Sep)' In Fig. ^ 
this means that one obtains the sensitivity reach by shift- 
ing the ellipse down on the pencil until the error bars of 
a fixed point (fixed #cpO touch the opposite-sign pen- 
cil. The sensitivity reach can then be read off as the 
corresponding value of sin 2 29^ for which this ellipse is 
drawn. Of course, this sensitivity reach depends on the 
true value of Sep, which determines the position on the 
ellipse. One can easily read off the bi-rate graph, that the 
best sensitivity reach is obtained close to Sep = 37r/2 and 
the worst close to Sep = tt/2. This behavior can also be 
found in the complete simulation result - it originates in 
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the strong dependence of the total event rates on the true 
value of Sep- In addition, the estimated sensitivities from 
the error bars in the bi-rate picture are a good estimate 
for the actual sensitivity reaches. For the negative-sign 
pencil, the ellipses are (for the same values of sin 2 2613) 
somewhat closer to the origin, which means that the sen- 
sitivity to the inverted mass hierarchy is somewhat worse. 
Nevertheless, the qualitative behavior as a function of the 
true value of Sep can be rediscovered. 

There are two more interesting aspects in Fig.^ First, 
since we do not know the true value of Sep realized by 
Nature, it is often convenient to show the conservative 
sensitivity reach in more condensed plots (c/., vertical 
lines in right panel). This conservative reach then de- 
scribes what region the experiment will be able to cover 
for sure, and helps to optimize it in terms of the other 
parameter values (such as in Figs. 4 and 8 of Ref. H^). 
Second, as we have discussed in Sec. |D] the size of Am|i 
is proportional to the sizes of the ellipses. Therefore, 
for a larger (smaller) Am^ within the solar- allowed re- 
gion, the ellipses become more (less) extended and the 
dependence on the true value of <5cp becomes stronger 
(smaller). Therefore, it can be shown that the amplitudes 
of the 5cp-dependencies in Fig.^ right panel, depend on 
Am^. However, mainly the conservative reaches (close 
to 3w/2) are affected for large Ato^ij whereas the best- 
case reaches (close to 7r/2) remain almost unaffected. 



V. THE SENSITIVITY TO CP VIOLATION 

We define sensitivity to CP violation if the experiment 
can distinguish Sep from the CP conserving values and 
tt at the chosen confidence level. This means that for a 
fixed simulated/true value of Sep ^ {0, tt}, the CP vio- 
lation sensitivity will be destroyed if any of the CP con- 
serving values or 7r fits the chosen simulated value <5cp ■ 
Since, in principle, any value 5cp ^ {Oj 71 "} violates CP, 
one can discuss the CP violation sensitivity as a function 
of this simulated value. Obviously, the closer Sep is to 
a CP conserving value or tt, the harder it will be to 
establish the (in this case small) CP violation. A special 
case of the sensitivity to CP violation is the sensitivity 
to maximal CP violation, which corresponds to the true 
values 5qp = tt/2 or 3tt/2. Since these values are farest 
away from CP conservation, the CP violation will in most 
cases be easiest to detect. In addition, we will refer to 
the CP violation for any simulated value 5cp 7^ 0, 7r as 
the sensitivity to any CP violation. 

Because of the correlation between sin 2 29i3 and <5cp , 
it is necessary to have some information on the CP con- 
jugated channels in order to efficiently measure Sap- It 
has, for example, been demonstrated in Refs. |29l l30l l32l ] 
that this implies a substantial fraction of antineutrino 
running in one or both of the JHF-SK and NuMI 
off- axis experiments. Alternatively, a large reactor ex- 
periment could s upp ly the CP-complementary informa- 
tion [l7l nil I42I l50| . However, it has been shown that 



CP violation measurements are very difficult for first- 
generation superbeams at the LMA-I solution. As we 
will demonstrate later, one of the best options for the 
superbeams is to operate JHF-SK with neutrinos and 
NuMI with antineutrinos only. Therefore, we use this 
configuration in this section as an example, and we will 
discuss other possibilities in the next section. In fact, we 
will observe that none of the other alternatives can do a 
CP violation measurements. 

We show in Fig.[21the potential to establish any CP vi- 
olation for JHF-SK (neutrino running mode) and NuMI 
(antineutrino running mode), i.e., the potential to estab- 
lish that a chosen value <$cp can be distiguished from 
the CP conserving solutions and tt. In the left-hand 
panel, the bi-rate graph is shown for the total charged 
current event rates of the two experiments, where the 
smaller antineutrino cross-section is compensated by the 
large NuMI detector in order to have similar error bars 
with similar statistical weights. In the right-hand panel, 
the region of sensitivity to CP violation (within thick 
black curve) is given as a function of the true values of 
sin 2 2#i3 and Sep for the complete statistical experiment 
simulation including all information. One can read off 
the sensitivity to maximal CP violation Sep = tt/2 or 
3tt/2 along the horizontal dashed lines in the right-hand 
plot. 

Let us first of all discuss what one could learn about 
the sensitivity to maximal CP violation from the bi- 
rate graph in Fig. [3 where we include the effect of the 
sgn(Am§ 1 )-degeneracy. For example, if we assume a nor- 
mal mass hierarchy and Sep = +tt/2 (stars), we will only 
be able to discover CP violation (at the chosen confidence 
level) if the error bars do not overlap the CP conserv- 
ing values and tt (triangles). Thus, the CP violation 
sensitivity is characterized by making the CP conserving 
values and tt "special" and it is determined by the rela- 
tive positions of the CP violating and CP conserving val- 
ues considered. It is often not directly representative for 
a CP precision measurement which does not emphasize 
any special values of Sep, and vice versa. For example, a 
CP precision of 180° at <5cp — tt/2 does not necessarily 
mean that one can establish maximal CP violation, since 
the fit manifold may not be exactly symmetric around 
Sep — tt/2 and one of the CP conserving solutions may 
be overlapped. 

The definition of the CP violation sensitivity also im- 
plies that overlapping one of the CP conserving values 
of any degenerate solution (which is present at the cho- 
sen confidence level) would destroy this sensitivity. For 
the sgn(Am| 1 )-degeneracy and a normal mass hierar- 
chy, this corresponds to an overlap with the triangles of 
the dashed curves in the bi-rate graph. In our example, 
the negative-sign ellipses are moved towards <5cp — tt/2 
and always from Sep — Stt/2 on the positive-sign el- 
lipses. This means that for ^cp = tt/2, the degenerate 
solution determines the CP violation sensitivity and for 
Sep = 3tt/2, the original solution (for the normal mass 
hierarchy). Thus, the degenerate solution is only im- 
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FIG. 2: (Color online) The sensitivity to any CP violation (<5cp 7^ 0, it) for JHF-SK (neutrino running only) and NuMI 
(antineutrino running only). The left-hand plot shows the bi-rate graph, where the ellipses are drawn for sin 2 20n = 0.03, 0.06, 
and 0.1, respectively, from the smallest to the largest. The right-hand plot marks the regions of sensitivity to CP violation as 
a function of the true values of sin 2 28n and 5cp for the complete statistical experiment simulation including all information, 
where the sensitivity is given within the gray-shaded region. The dashed curve refers to the additional sensitivity if one is not 
taking into account the degenerate solution. The vertical lines correspond to the values of sin 2 2#i3 for which the ellipses in the 
left-hand plot are drawn. The horizontal lines are marked by the corresponding symbols from the left-hand plot, where solid 
lines refer to CP conservation and shaded lines to maximal CP violation. The symbols and other parameters are the same as 
in Fig. but we assume a normal mass hierarchy in the right-hand plot. For the inverted mass hierarchy, the role of the two 
peaks in the right-hand plot would be exchanged. 



portant for the sensitivity to Sep = 7r/2. In addition, 
this interplay between the original and degenerate solu- 
tions can easily be read off the bi-rate graphs for the in- 
verted mass hierarchy and other setups. For example, for 
the inverted mass hierarchy, the role of <5cp — tt/2 and 
5qp — 3tt/2 is just exchanged. However, it should be 
noted that for high-luminosity upgrades, the degenerate 
solution does not necessarily appear at all below the cho- 
sen confidence level. For the first-generation superbeams, 
it is normally present at very low x 2 -values. Eventually, 
one can immediately read off the bi-rate graph, that the 
sensitivity to maximal CP violation is given for large val- 
ues for sin 2 2#i3 close to Sep = 3n/2 only, which is con- 
sistent with the complete analysis in the right-hand panel 
of Fig. |21 (along horizontal dashed lines). 

One can also estimate the sensitivity to any CP vio- 
lation Sep 7^ 0, 7r by following the points on a specific 
ellipse in the bi-rate picture until the error bars touch 
one of the CP conserving solutions (triangles) . With this 
approach, one can read off how close one can get to the 
CP conserving solutions. Indeed, the estimates from the 
bi-rate picture are consistent with the result from the 
complete experiment simulation in the right-hand panel 
of Fig. |21 One can easily see by the comparison of the 
dashed and thick solid curves in the right-hand panel that 
the impact of the sgn(Arn| 1 )-degeneracy is very large, 



since the strong matter effects at the NuMI beam cause 
a large separation of the positive- and negative-sign pen- 
cils. It completely destroys the sensitivity to maximal 
CP violation Sep = tt/2, while the sensitivity to maximal 
CP violation 3tt/2 is still given down to sin 2 2#i3 ~ 0.03. 
In addition, the sensitivity to any CP violation is only 
present very close to the maximally CP violating value. 
Again, one can read off this figure that the sgn(Am| 1 )- 
degeneracy only affects the first and second quadrants 
(for the positive mass hierarchy) |2^| , since the negative- 
sign ellipses are shifted into the direction of <5cp = tt/2 
for the considered setup. Thus, the CP conserving val- 
ues are shifted towards Sep = tt/2 and make the sensi- 
tivity to CP violation worse for any value within these 
quadrants. From the dashed curve in the right-hand 
figure without sgn(Arn| 1 )-degeneracy, it is obvious that 
the sensitivity to any CP violation is symmetric before 
the inclusion of this degeneracy. Thus, in this case, the 
sgn(Am§ 1 )-degeneracy is the main reason for the asym- 
metry of the CP sensitivities between the first/second 
and third/fourth quadrants. 

The concepts in Fig. can also be extended to super- 
beam upgrades, such as JHF-HK. In fact, one can also 
observe that only the first and second quadrants (around 
Sqp = tt/2) are affected by the sgn(Am| 1 )-degeneracy in 
a similar way to this figure. For JHF-HK, correspond- 
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ing figures can be found as Fig. 18 in Ref. [2lJ or Fig. 7 
in Ref. |&l|. Note that for neutrino factories, the behav- 
ior is somewhat different and more complicated, because 
they use a much broader energy spectrum and a different 
L/E. 



VI. THE PRECISION OF <5 CP : 
"<5cp-PATTERNS" 

We define the precision of Sep as the fraction of the 
CP circle which fits the selected true value of Sep at the 
chosen confidence level. Therefore, we also call it "CP 
coverage" , since it describes what fraction of the fit val- 
ues is "covered" on the CP circle by an individual mea- 
surement, and allows to include complicated topologies 
in a very simple way. A coverage (or precision) of 360° 
corresponds to all values of 5cp fitting the true value (and 
thus, no information on 5qp) and a coverage (or preci- 
sion) close to 0° to a very precise measurement. Hence, 
the value of 360° minus the CP coverage corresponds to 
the range which can be excluded. Any value of the preci- 
sion which is smaller than 360° will therefore imply some 
information on i5cp and exclude certain regions on the 
CP circle. From this definition, the difference to the CP 
violation sensitivities is obvious: The CP precision mea- 
surement answers the question if we can learn something 
about 5cp at all, and the CP violation measurement if 
we can distinguish CP violation from CP conservation. 
For example, for the case of Sep close to or ir, there is 
no sensitivity to CP violation, but we may still be able 
to exclude certain regions on the CP circle. For the CP 
precision measurement, the treatment of the sgn(Amg 1 )- 
degeneracy is straightforward and is given by its defini- 
tion: Any solution of the opposite-sign circle fitting the 
true value of Sep at the chosen confidence level has to be 
removed from the allowed regions on the CP circle. 

In Fig. |3| we show the CP coverage as a function of 
the true value of <5cp for a complete statistical analy- 
sis including all information and different combinations 
of next-generation experiments as examples. Obviously, 
the precision strongly depends on the true value of Sep • 
Especially the dependence on the true value of Sep allows 
to quote a wide range of values for the CP precision of a 
specific experiment or combination of experiments. One 
can also see that this dependence is very complicated, 
although there are some regularities. We further on refer 
to these regularities as a function of the true value of Sep 
as "CP patterns" . 

Let us first compare the CP coverage measurement in 
Fig. El to the CP violation measurement. In order to 
have sensitivity to maximal CP violation, it is necessary 
to have a CP precision smaller than 180° at the CP vi- 
olating values Sep = n/2 or 3tt/2. 2 Thus, only the first 



2 It is necessary, but not sufficient, since the fit manifold may not 



experiment combination (JHF-SK with neutrinos and 
NuMI with the large detector and antineutrinos) could 
have some sensitivity to maximal CP violation (indeed, 
it has, as it can be seen from Fig. EJ. All the other com- 
binations do not have a sensitivity to CP violation, as 
it has been demonstrated in Ref. j2(| for combinations 
of first-generation superbeams. However, Fig. [31 shows 
that one may still learn something about the CP phase 
if sin 2 2#i3 is large, which strongly depends on the true 
value of Sep- In order to exclude certain ranges for (5cp, 
the CP coverage has to be smaller than 360°. Since we 
do not know yet what the true value of <5cp is, one can 
be lucky or unlucky if one is building an experiment with 
large ranges with a CP coverage of 360°. Thus, the risk 
for such an experiment not to supply any new informa- 
tion is rather high. 

In the following, we will give a qualitative description 
of the CP patterns of the next-generation experiments by 
comparing Fig. El upper left panel, with Fig. El left, since 
these two figures are computed for the same experiments 
and parameter values. Note that the largest ellipses in 
Fig. computed for sin 2 28i 3 = 0.1 correspond to Fig. El 
The other experiments in Fig. [31 show the same qualita- 
tive behavior. 3 In order to keep this description simple, it 
is first of all useful to ignore the effect of the sgn(Am§ 1 )- 
degeneracy, i.e., to concentrate on the dashed curves. In 
Fig. El left, the error bars are relatively large compared to 
the size of the ellipses. In addition, systematics and cor- 
relations make them even larger than they appear there. 
Thus, for Sep = or tt the error bars cover the whole el- 
lipse and one cannot obtain any information on i5cp, i.e., 
the error in Fig. El (left) reaches 360°. In particular, they 
always touch the opposite side of the ellipse, because the 
correlation with sin 2 2#i3 allows an ellipse for a different 
value of sin 2 2$i3 to overlap the original one. However, 
for Sep = tt/2 or 37r/2 only half of the error bars overlap 
the ellipse, which means that there the information on 
Sep is significantly better. Depending on the size of the 
error bars, and thus the true value of sin 2 26*i3, the CP 
coverage can differ up to a factor of two between Sep = 
or 7T and Sep = 7r/2 or 37r/2. 

As far as the sgn(Am§ 1 )-degeneracy is concerned, we 
illustrate its effects for a normal mass hierarchy by the 
dashed (no degeneracy) and solid (final result) curves in 
Fig. El Since, in this example, the error bars are rather 
large, it is only effective for the first and second quadrants 
in a straightforward way, because there "new" values of 
Sep enter the error bars by the shifting of the negative 
mass hierarchy ellipse (c/., Fig. El left). Because the pos- 



be symmetric around the maximally CP violating value and still 
overlap one of the CP conserving solutions. 
3 Especially, similar to the CP-conjugated channel at a superbeam, 
the reactor experiment can also help to resolve the correlation 
between sin 2 20i3 and Sep with a precise sin - ' 29m determina- 
tion ll ^ , which means that one would not expect large surprises 
in this case. 
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FIG. 3: The CP coverage as a function of the true value of <5cp for the four different experiment combinations as given in the 
plot labels and a complete statistical experiment simulation including all information (90% confidence level): JHF-SK with 
neutrino running plus NuMI with antineutrino running, JHF-SK CC in the combined neutrino-antineutrino running mode (25% 
of the total running time in the neutrino mode and 75% in the antineutrino running mode), JHF-SK plus Reactor-II, and 
NuMI plus Reactor-II. Here Reactor-II refers to a large reactor experiment with near and far detectors, as it is defined 
in Ref. |lq| with an integrated luminosity of 8 000tGWy. In these plots, the CP coverage [Degrees] is defined as the part of 
the CP circle which fits the true value which is given on the horizontal axis, i.e., 360° corresponds to no information on 5cp- 
For the oscillation parameters the current best-fit values with sin 2 2#i3 = 0.1 and a normal mass hierarchy are assumed. The 
dashed curves refer to not taking into account degeneracies. 



itive and negative mass hierarchy ellipses are rather dif- 
ferent from each other for large values of sin 2 2#i3, but 
quite similar to each other for small values of sin 2 2#i3, 
the effects of the degeneracy are much larger for the large 
values of sin 2 2#i3 (provided that it is present at all at the 
chosen confidence level). For the inverted mass hierar- 
chy, one would just observe the inverted behavior for the 
exchanged quadrants. Thus, the sgn(Am§ 1 )-degeneracy 
breaks the symmetry between 5cp = n/2 and 37r/2 [23l| . 



Compared to exclusion measurements of the leptonic 
CP phase, it is interesting to investigate high-precision 
measurements, such as from a superbeam upgrade. We 
here focus on JHF-HK, which has a very narrow energy 
spectrum and is therefore well-suited for the bi-rate ap- 
proach. In addition, the CP coverage is for very precise 
measurements identical to the precision and can there- 
fore describe both exclusion and high-precision measure- 
ments. Fig. 0] shows the CP precision as a function of 
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FIG. 4: The CP coverage (upper row) as a function of the true value of Sep for JHF-HK and three different values of sin 2 2613 
as given in the plot labels (90% confidence level). The plots in the upper row include the complete statistical simulation 
including all information and correspond to the ellipses in the lower row. For the oscillation parameters the current best-fit 
values and a normal mass hierarchy are assumed. The dashed curves refer to not taking into account degeneracies. The meaning 
of the symbols is explained in Fig. 



<5cp (upper row) for three different values of sin 2 26*13 (in 
columns). In the lower row of Fig. the corresponding 
ellipses are drawn. Note that JHF-HK is operated in the 
combined neutrino-antineutrino running mode with two 
years of neutrino running and six years of antincutrino 
running. 

For JHF-HK, the statistical error bars are very small 
compared to the ellipses. Thus, without taking into ac- 
count degeneracies (dashed curves in upper row), the CP 
patterns are, in principle, very regular (though differ- 
ent from the case of large error bars - see below) . How- 
ever, there are some complications. First, for sin 2 26*13 = 
0.001, the superbeam upgrade is dominated by system- 
atics, which means that the actual error bars (including 
both statistics and systematics) are much larger than the 
statistical error bars in the figure. Second, the effects of 
the degenerate solution can either be rather large and 
may lead to very irregular patterns (e.g., for sin 2 26*13 = 
0.01) or the degeneracy may not be present at all at the 
chosen confidence level (e.g., for sin 2 26*13 = 0.1). Com- 
paring the ellipses to the corresponding CP pattern, we 
can distinguish three cases with respect to the size of the 
error bars and the sgn(Am| 1 )-degeneracy: 



Very small values sin 2 2#i3. For small values of 
sin 2 26*13 (left-hand panel of Fig. |3J| , the actual 
error bars are much larger that the statistical error 
bars in Fig. 0] (left) because of the systematics 
domination of the background and its uncertainty, 
as well as other effects. Thus, we have a similar 
case to the first-generation superbeams, where 
the CP coverage is best close to 7r/2 and 3/27T 
(without degeneracy). In addition, the sgn(Am| 1 )- 
degeneracy is present at the chosen confidence 
level. However, the ellipses for the original and 
degenerate solutions are very similar to each 
other and the actual error bars are rather large, 
which means that the effect of the degeneracy is 
small. Note that this specific case is an example 
of the limitations of the bi-rate approach, since 
the systematics is more complicated to be simply 
included in the bi-rate picture. 

Large values of sin 2 2#i3. For large values of sin 2 26 ) i3 
(see right column of Fig. |3J|, i.e., sin 2 26>i 3 > 0.01, 
we obtain the simplest pattern with a continuous 
dependence on Sqp- Since in this case the error 
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bars are very small compared to the size of the el- 
lipses, the CP precision becomes best at and w 
and worst at ir/2 and 3/2ir because of the shallower 
dependence on Sep close to and ir. In addition, 
the sgn(Am§ 1 )-degeneracy is normally not present 
at the chosen confidence level. Thus, it hardly af- 
fects the CP pattern. 

Intermediate values of sin 2 29±3. For intermediate 
values of sin 2 20i3 (see middle column of Fig. 0}, 
the pattern without degeneracy correspond to 
the one for large sin 2 20i3 and is very regular 
(c/., dashed curve). However, the degenerate 
solution is present at the chosen confidence level 
and can have a very large effect especiall y fo r 

0. 001 < sin 2 20i3 < 0.01 (c/., Fig. 8 of Ref. :2lj). 
Although its effect is, in principle, stronger in 
the first and second quadrants around Sep = vr/2 
than in the third and fourth quadrants (such 
as for JHF-SK), it can also affect the other 
quadrants because the ellipses for the original 
and degenerate solutions are still very close to 
each other in this sin 2 20i3-range. Since the error 
bars are very small compared to the size of the 
ellipse, the CP coverage can be easily doubled if 
the degenerate solution moves "new" values of 
Sqp into the error bars of the original solution, 

1. e., the CP precision becomes worse. However, 
in some cases (shortly after and close to 3ir/2) 
the values of the negative-sign ellipse correspond 
to the ones of the positive-sign ellipse and the 
degenerate solution has no effect. 4 In other cases, 
only a fraction of new values is moved into the 
error bars of the original solution, which means 
that the CP coverage is in-between the dashed 
curves (in the upper row of Fig. [3} and their 
doubles. This complex interaction of the original 
and opposite-sign solution leads to the irregular 
CP pattern after inclusion of the degeneracy for 
sin 2 26>i3 = 0.01. 

In summary, the CP patterns for JHF-HK are much 
more complicated than the ones for the first-generation 
superbeams. However, resolving the sgn(Am§ 1 )- 
degeneracy would in this case help the CP performance 
very much and could even double the precision for inter- 
mediate values of sin 2 29^. 



VII. SUMMARY AND CONCLUSIONS 

In this study, we have investigated the mass hierarchy, 
CP violation, and CP precision measurements at neu- 



4 Note that again the opposite side of the ellipse with respect to 
sin 2 2613 has to be included because of the correlation between 
5cp and sin 2 29\:j. 



trino superbeams as function of the true value of Sep- 
We have used a likelihood analysis including system- 
atics, correlations, degeneracies, spectral information, 
and a realistic detector simulation. Since the results 
from such an analysis are very difficult to interpret, we 
have chosen bi-rate graphs to illustrate many CP-phase 
dependent properties. Bi-rate graphs are based upon 
bi-probability graphs, but they use the total charged- 
current appearance rates instead of the oscillation prob- 
abilities in the bi-probability graphs. Since error bars 
can be attached to each point in a bi-rate graph, one can 
directly read off quantitative estimates for many effects. 
Furthermore, they are especially useful for off-axis super- 
beams, because these are "quasi-monochromatic" narrow 
band beams, which only carry little spectral information. 
Therefore, to a first approximation, only one ellipse in the 
bi-rate picture corresponds to one choice of energy and 
baseline. By comparing the bi-rate graphs with the com- 
plete statistical simulation, we have demonstrated that 
the role of the true value of Sep in the mass hierarchy, CP 
violation, and CP precision measurements at superbeams 
can be understood in terms of these graphs. Though it 
can be used as a powerful tool in many cases, the bi-rate 
picture does not always exactly fit the complete calcu- 
lations, since in some cases there are significant contri- 
butions of some spectral information, the disappearance 
channels, systematics, and correlations other than be- 
tween sin 2 20i3 and Sep ■ We have shown these limita- 
tions at the example of the systematics domination of 
JHF-HK for small values of sin 2 20 i3 . 

For the mass hierarchy determination, we have demon- 
strated that the true value of <5cp affects the potential to 
measure the mass hierarchy in a way that the sin 2 20i3- 
sensitivity reach is worst close to Sep — 7r/2 (conser- 
vative reach) and best close to 5cp = 37r/2. Further- 
more, the best case sin 2 20i 3 -sensitivity reach (close to 
Sep = 37r/2) hardly depends on the true value of Ato 2 i, 
whereas the conservative case sin 2 20i3-sensitivity reach 
(close to Sep = tt/2) is highly affected by large values of 
Amli- Since the true value of <5cp is given by Nature, 
it is often convenient to use the conservative reach for 
Sep — tt/2 in order to determine what an experiment can 
do in either case. Eventually, the sensitivity to the nor- 
mal mass hierarchy is better than the one for the inverted 
mass hierarchy. All of these effects can be qualitatively 
demonstrated with the help of bi-rate graphs. 

The sensitivity to CP violation depends on the relative 
positions of the original and degenerate CP conserving 
values (0 and 7r) compared with the considered CP vio- 
lating value. If the error bars at a CP violating solution 
in the bi-rate graph overlap one of the CP conserving 
solutions, then the sensitivity to CP violation will be de- 
stroyed. The asymmetry of the CP violation sensitivity 
between the first/second Sep G]0, 7t[ and third/fourth 
quadrants <5cp &]tt,2tt[ is caused by the sgn(Am§i)- 
degeneracy. For the normal mass hierarchy, primarily 
the sensitivities to CP violation in Sep €]0,7r[ are af- 
fected by the sgn(Am 2 1 )-degeneracy, and for the inverted 
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mass hierarchy, primarily the sensitivities to CP violation 
in Sqp €]tt, 2tt[. Note that CP violation measurements 
are very difficult for first-generation superbeams. In fact, 
only the combination of JHF-SK (neutrino running only) 
combined with NuMI (antineutrino running only) would 
only have a sensitivity close to maximal CP violation 
<5cp — 7r/2 (for inverted mass hierarchy) or <5cp — 3ir/2 
(for normal mass hierarchy). 

Unlike the CP violation measurement, the CP preci- 
sion measurement indicates if one could learn something 
about Sep at all. For example, one may not be able to 
detect CP violation, since Sep is too close to or w, but 
one may still be able to exclude certain ranges for Sep- 
We have used the "CP coverage" as a useful quantity 
to describe the 5cp-performance of an experiment. It 
is the range of ^cp-values (as fraction of the CP circle) 
which fit a chosen simulated value (true value) of Sep- 
Therefore, the concept of the CP coverage takes into ac- 
count the cyclicity of the Jcp-dependence. In addition, 
it can describe <5cp -exclusion measurements, as well as 
high-precision measurements. The CP coverage strongly 
depends on the true value of 5cp itself, and this depen- 
dency results in rather complicated "CP patterns" . The 
CP patterns depend on sin 2 2#i3 and the presence of the 
sgn(Amg 1 )-dcgeneracy, and can be qualitatively under- 
stood in terms of bi-rate graphs. The strong dependence 
on the true values of Sep itself implies that one has to 
be careful when quoting or comparing CP precision per- 
formances of different experiments. For example, an ex- 
periment, such as a first-generation superbeam, may not 
supply any new information on <5cp in a wide range of 
true values of <5cp- Thus, it is always useful to show and 
compare the complete CP patterns for different exper- 
iments. In particular, we have also demonstrated that 
for the normal mass hierarchy, only the precisions for 
Sqp €]0, 7t[ are affected by the sgn(Am| 1 )-degeneracy, 
and for the inverted mass hierarchy, only the precisions 
for Sqp € ]?r, 2*7r[ for the first-generation superbeams. 
Thus, the sg^AmJ-jJ-degeneracy causes an asymmetry 
between the first/second Sep €]0, n[ and third/fourth 
quadrants <5cp € ]n, 2ir[. For the superbeam upgrades, 
the situation is more complicated, since the error bars 
become very small and the sgn(Am§ 1 )-degeneracy may 



or may not be present at the chosen confidence level. 
Especially for intermediate values of sin 2 26*i3, the pres- 
ence of the sgn(Am 2 1 )-degeneracy causes very irregular 
CP patterns. In addition, especially for high precision 
instruments, the appearance of all degeneracies as func- 
tions of the confidence level is the reason that the CP 
precision over-proportionally depends on the chosen con- 
fidence level. Thus, it is important that the same confi- 
dence level be used for a comparison of different experi- 
ments. 

Finally, bi-rate graphs have been demonstrated to be a 
useful tool, which links the qualitative predictions of the 
bi-probability picture with complete quantitative evalua- 
tions including spectral information, systematics, corre- 
lations, and degeneracies. For example, the complicated 
CP patterns of the CP precision measurements can be un- 
derstood in terms of bi-rate graphs. Of course, it would 
also be useful to extend this discussion to neutrino facto- 
ries. However, neutrino factories have a lot more spectral 
information than superbeams, and therefore, their treat- 
ment in terms of bi-rate diagrams will be much more 
complicated. Nevertheless, one can still find many of the 
qualitative features observed in this work in neutrino fac- 
tory measurements. 



Acknowledgments 

I would like to thank Steve Geer, Deborah Harris, 
Patrick Huber, Manfred Lindner, and Tommy Ohlsson 
for useful discussions. In addition, I especially would 
like to thank Stephen Parke for extensive discussions and 
valuable comments, and the Fermilab theory group for 
their warm hospitality during my stay within the sum- 
mer visitor program. 

This work has been supported by the "Sonder- 
forschungsbereich 375 fur Astro- Tcilchcnphysik dcr 
Deutschen Forschungsgemeinschaft" (SFB 375 of DFG), 
the "Studienstiftung des deutschen Volkes" (German Na- 
tional Academic Foundation), and the Fermilab summer 
visitor program. 



[1] Y. Fukuda et al. (Supcr-Kamiokande Collaboration), 

Phys. Rev. Lett. 81, 1562 (1998), hep-ex/9807003. 
[2] Y. Fukuda et al. (Super-Kamiokande Collaboration), 

Phys. Rev. Lett. 82, 2644 (1999), hep-ex/9812014. 
[3] T. Toshito (Super-Kamiokande) hep-ex/0105023. 
[4] Q. R. Ahmad et al. (SNO Collaboration), Phys. Rev. 

Lett. 87, 071301 (2001), nucl-ex/0106015. 
[5] Q. R. Ahmad et al. (SNO Collaboration), Phys. Rev. 

Lett. 89, 011301 (2002), nucl-ex/0204008. 
[6] S. N. Ahmed et al. (SNO) nucl-ex/0309004. 
[7] K. Eguchi et al. (KamLAND Collaboration), Phys. Rev. 

Lett. 90, 021802 (2003), hep-ex/0212021. 



[8] M. Apollonio et al. (CHOOZ Collaboration), Phys. Lett. 

B466, 415 (1999), hep-ex/9907037. 
[9] M. Apollonio et al. (CHOOZ Collaboration), Eur. Phys. 

J. C27, 331 (2003), hep-ex/0301017. 
[10] K. Nakamura (K2K Collaboration), Nucl. Phys. Proc. 

Suppl. 91, 203 (2001). 
[11] V. Paolone (MINOS Collaboration), Nucl. Phys. Proc. 

Suppl. 100, 197 (2001). 
[12] D. Duchesneau (OPERA), eConf C0209101, TH09 

(2002), hep-ex/0209082. 
[13] Y. Itow et al, Nucl. Phys. Proc. Suppl. Ill, 146 (2001), 

hep-ex/0106019. 



14 



[14] D. Ayres et al. hep-ex/0210005. 
[15] D. Beavis et al. hep-ex/0205040. 

[16] V. Martemyanov, L. Mikaelyan, V. Sinev, V. Kopeikin, 

and Y. Kozlov hep-ex/0211070. 
[17] H. Minakata, H. Sugiyama, O. Yasuda, K. Inoue, and 

F. Suekane, Phys. Rev. D68, 033017 (2003), hep- 

ph/0211111. 

[18] P. Huber, M. Lindner, T. Schwetz, and W. Winter, Nucl. 

Phys. B665, 487 (2003), hep-ph/0303232. 
[19] M. H. Shaevitz and J. M. Link hep-ex/0306031. 
[20] M. Apollonio et al. and references therein, hep- 

ph/0210192. 

[21] P. Huber, M. Lindner, and W. Winter, Nucl. Phys. 

B645, 3 (2002), hep-ph/0204352. 
[22] J. Burguet-Castell, M. B. Gavela, J. J. Gomez-Cadenas, 

P. Hernandez, and O. Mena, Nucl. Phys. B608, 301 

(2001) , hep-ph/0103258. 

[23] H. Minakata and H. Nunokawa, JHEP 10, 001 (2001), 

hep-ph/0108085. 
[24] G. L. Fogli and E. Lisi, Phys. Rev. D54, 3667 (1996), 

hep-ph/9604415. 
[25] V. Barger, D. Marfatia, and K. Whisnant, Phys. Rev. 

D65 , 073023 (2002), hep-ph/0112119. 
[26] Y. F. Wang, K. Whisnant, Z.-h. Xiong, J. M. Yang, and 

B.-L. Young (VLBL Study Group H2B-4), Phys. Rev. 

D65, 073021 (2002), hep-ph/0111317. 
[27] K. Whisnant, J. M. Yang, and B.-L. Young, Phys. Rev. 

D67, 013004 (2003), hep-ph/0208193. 
[28] V. Barger, D. Marfatia, and K. Whisnant, Phys. Rev. 

D66, 053007 (2002), hep-ph/0206038. 
[29] P. Huber, M. Lindner, and W. Winter, Nucl. Phys. 

B654, 3 (2003), hep-ph/0211300. 
[30] H. Minakata, H. Nunokawa, and S. Parke, Phys. Rev. 

D68, 013010 (2003), hep-ph/0301210. 
[31] P. Migliozzi and F. Terranova, Phys. Lett. B563, 73 

(2003), hep-ph/0302274. 
[32] V. Barger, D. Marfatia, and K. Whisnant, Phys. Lett. 

B560, 75 (2003), hep-ph/02 10428. 
[33] J. Burguet-Castell, M. B. Gavela, J. J. Gomez-Cadenas, 

P. Hernandez, and O. Mena, Nucl. Phys. B646, 301 

(2002) , hep-ph/0207080. 

[34] A. Donini, D. Meloni, and P. Migliozzi, Nucl. Phys. 

B646, 321 (2002), hep-ph/0206034. 
[35] D. Autiero et al. hep-ph/0305185. 
[36] A. Asratyan et al. hep-ex/0303023. 
[37] W. Winter hep-ph/0308227. 

[38] H. Minakata and H. Nunokawa, Nucl. Instrum. Meth. 

A503, 218 (2001), hep-ph/0111130. 
[39] H. Minakata, H. Nunokawa, and S. Parke, Phys. Lett. 

B537, 249 (2002), hep-ph/0204171. 
[40] H. Minakata, H. Nunokawa, and S. Parke, Phys. Rev. 

D66, 093012 (2002), hcp-ph/0208163. 
[41] H. Minakata, H. Nunokawa, and S. Parke, AIP Conf. 

Proc. 670, 132 (2003), hep-ph/0306221. 
[42] H. Minakata, H. Nunokawa, and S. Parke hep- 

ph/0310023. 

[43] A. Cervera et al, Nucl. Phys. B579, 17 (2000), hep- 
ph/0002108. 

[44] M. Freund, P. Huber, and M. Lindner, Nucl. Phys. B585, 

105 (2000), hep-ph/0004085. 
[45] M. Freund, Phys. Rev. D64, 053003 (2001), hep- 

ph/0103300. 

[46] V. D. Barger, D. Marfatia, and B. P. Wood, Phys. Lett. 
B498, 53 (2001), hep-ph/0011251. 



[47] M. C. Gonzalez-Garcia and C. Pena-Garay, Phys. Lett. 

B527, 199 (2002), hep-ph/0111432. 
[48] M. Maltoni, T. Schwetz, and J. W. F. Valle, Phys. Rev. 

D67, 093003 (2003), hep-ph/0212129. 
[49] M. Maltoni, T. Schwetz, M. A. Tortola, and J. W. F. 

Valle hep-ph/0309130. 
[50] H. Minakata and H. Sugiyama hep-ph/0309323. 
[51] T. Ohlsson and W. Winter, Phys. Rev. D68, 073007 

(2003), hep-ph/0307178. 



